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Abstract 

We argue that an oscillatory motion of impurity-pinned solitons may occur in layered antifer- 
romagnetic compounds. The characteristic frequencies of these modes, that may be detected by 
resonance or inelastic neutron scattering, are estimated analytically and depend on the soliton sizes 
and types. 
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In recent years there has been much interest in the static and dynamic properties of 
layered compounds. These materials are well-known to exhibit quasi- two-dimensional (2D) 
magnetism at low temperatures and, in addition to the possibilities in technological ap- 
plications, they have been useful in testing theories pertinent to low dimensional systems. 
The 2D behavior of these structures arises because the ratio of the interlayer coupling J z 
to intralayer coupling J is of the order J z /J ~ 10~ 3 — 10~ 5 . The 2D effects can be exper- 
imentally observed in a narrow temperature interval just above the ordering temperature. 
The theory of these systems is contained in the Heisenberg Hamiltonians (isotropic and 
anisotropic). The isotropic Heisenberg magnets in two dimensions are described by the 
Hamiltonian H = iJJ^ j ' ^j, where 5, is the spin located at site i, the summation is 
taken over the nearest-neighbor sites on a square lattice, and the - and + signals represent 
ferromagnetic and antiferromagnetic coupling constants respectively. The continuum ver- 
sion of this Hamiltonian is given by the 2D 0(3) non-linear a model, which has long been of 
interest also to particle theorists because of its strong analogy with four- dimensional Yang- 
Mills theory. For Heisenberg antiferromagnet (AFM), the non-linear a model is described 
in terms of the Neel state and gives the antiferromagnetic long-wavelength fluctuations in 
the limit of zero temperature. It is represented by the following Hamiltonian : 

H, = \j J [(d,n) ■ {d»n)]d 2 x, // = 1, 2, (1) 

where n denotes the Neel unit vector {n\ + n 2 . + n\ = 1) characterizing the direction of 
the sublattice magnetization of the antiferromagnet. It is also well-known that this model 
supports interesting soliton-like solutions, obtained by Belavin and Polyakov[I[] using topo- 
logical considerations. In terms of the conformal representation w = (ni + m 2 )/(l — n 3 ), 
which is simply the stereographic projection of the Neel vector onto the complex plane, the 
Belavin-Polyakov soliton (5F-soliton) located at the origin and with unitary topological 
charge is expressed as w s = z/R, where z = x + iy and R is the soliton size. It implies that 
the third Neel vector component n% is given by 

Furthermore, it is useful to change variables from (x,y) to {n%,<f> = tan _1 (y/x)), in terms 
of which we can express the ratio r^/ni = tan0. Qualitatively, the configuration is radial, 
pointing down at the origin and up at infinity. The energy of these excitations is found to 



be independent of the soliton size due to the scale invariance of the non-linear a model and 
it is given by E s = 4irJ. The significance of these non-linear excitations for 2D magnets was 
recognized early in connection with the critical properties of these systems. For example, it 
was shown that the existence of large localized excitations will cause the correlation length 
to remain finite at any nonzero temperature JTJ . Considering the dynamical properties, the 
translation motion of this structure may be responsible for a central peak in the out-of-plane 
dynamic correlation function^, [|. Recently, Zaspel et a/@, ||| have shown that the Belavin- 
Polyakov solitons dominate the thermodynamics in the fluctuation region immediately above 
the Neel temperature of a large class of nearly classical 2D antiferromagnets. Besides, these 
authors || [7j have also shown that the introduction of a very small amount of nonmagnetic 
impurities into the magnetic sites of a classical 2D antiferromagnet creates a new type of 
static (impurity-pinned) soliton that affects the Arrhenius, e~ E ^ T ', temperature-dependent 
EPR linewidth by drastically changing the parameter E s . Using an approach on a discrete 
lattice, they considered the spin vacancy located at the soliton center and obtained two 
different impurity-solitons ( the P-soliton and the /-soliton), which were detected in the 
layered antiferromagnet {C^HjN H^)2M x Mni- x Cl4 : through EPR measurements |6|, 0]. The 
.P-soliton has the same structure of the PF-soliton but without a spin at its center, which 
results in a vortex singularity at the center. The /-soliton has smaller size and energy than 
the P-soliton, so it also has a configuration slightly different from the P-soliton. In this letter 
we argue that these solitons oscillate on the static spin vacancy and estimate the frequencies 
of these oscillatory modes of solitons analytically. It is suggested that this effect may be 
observed in electron spin resonance or inelastic neutron scattering experiments. 

In a recent paper, Pereira and Pires|| considered the problem of a soliton interacting with 
a nonmagnetic impurity in the continuum limit by modifying Hamiltonian (1) introducing 
a nonmagnetic impurity potential as follows: 

Hi = \j J [(d,n) ■ (d»n)}V(r)d 2 x, fi = l,2, (3) 

where V(f) = 1 if \r — d\ > a, and V(r) = if \r — d\ < a. Here, the impurity is localized at 
position d from the origin, and a is the lattice constant. This lack of magnetic interaction 
inside the circle of radius a, means that a spin located at d was removed from the lattice. 
From now we will consider the nonmagnetic impurity localized at the origin, d = (0, 0). If 
the soliton center coincides with the impurity center, this field theory reproduces the two 
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soliton solutions of Refs. || and still gives the sizes, energies and configurations of these 
solitons in an analytical way. For the P-soliton, almost entire space is mapped onto the 
plane, but part of the spin space sphere around the origin is not included in the mapping. It 
means that the stereographic projection of the Neel vector onto the complex plane is perfect 
for r = \z\ > a, so we can write the P-soliton configuration as 

n 2P J for r < a, 

= tan0; n 3P = < (4) 

nip I 77-3 for r > a. 

For the /-soliton, the mapping is smoothly deformed, leading to the following configuration 

n 2/ J for r < a, 

= tan</>; n 3I = < (5) 

n n I n 3 + &(P)(l -n|) for r > a, 

where the function b(R) was obtained minimizing the energy of excitation (5) for each value 
of the soliton size and it is given in the interval < b(R) < 1/2||. The subscript P and 
/ refer to the topological and impurity-centered solitons (P- and /-solitons) respectively. 
The energy of these two solitons depends on R with the /-soliton energy smaller than the 
P-soliton energy, but as R — > oo, Ej — > Ep — > 47r J. This continuum theory permits also to 
compare the solitons size (P/ = [(1 + 4b 2 ) 1 / 2 — 2b]R, Rp = R) showing that the /-soliton 
is smaller than the P-soliton||. 

The field theory given by (3) still shows that, if the soliton center is displaced in relation to 
the nonmagnetic impurity center, an effective interaction between the soliton and impurity 
that is attractive at short range is induced. If the soliton center and the spin vacancy center 
are separated by a distance r , the system energy can be approximated by|8[ 

2 

4 \r 2 + R 2 

where 



J - J ro>a — - C/ s 



i-TlV^l +P-(r ,P)-P + (r ,P) 



(6) 



1 ' ' [(r ± a) 2 + r 2 ]V2[( ro ± G ) 4 - P 4 ] ' { } 

Equation (6) shows that, for r < P, the soliton energy decreases as the soliton center 
approximates to the impurity. It indicates that there is a short-range attractive well with 
length scale on the order of P. On the other hand, if the soliton center is localized at 
distance r > P from the spin vacancy, the effective potential between the soliton and the 
nonmagnetic impurity is repulsive. However, although expression (6) correctly suggests that 



d 2 E\ T ^ E S -EJR) 

= K = 2 ^ ^-L + 2E S 



(9) 



it is energetically favorable for the soliton to nucleate about the spin vacancy, this equation 
is valid only for separations larger than one lattice spacing (ro > In fact, in the limit 

that the separation vanishes (r — > 0), the soliton configuration changes to either the P- 
soliton or /-soliton configuration, which are radially symmetric and represent local minima 
of energy 0. Then, if there is a spin vacancy exactly at the soliton center, Eq.(6) should give 
the values Ej(R) or Ep(R), depending on which soliton was formed (the I- or P-soliton). 
Hence, in the limit ro — > 0, we use an expansion of the energy around the bottom of the well 
as follows 

E^in, «) S EaiR) + \rl (g)^ + ^ (f )^ + , (8, 

where a =P or / represents one of the two possible solitons. For small displacements, we 

can neglect the high-order terms in the expansion. The complete function E(ro, R) is likely 

to be a smooth interpolation between the limit (8) and Eq. (6), with the transition between 

these expressions occurring at r ~ a. Keeping only the harmonic term in (8), we get 

R 2 R 2 

_V5[R 4 - (2a) 4 ] 4(i? 2 + a 2 ) 2 _ 

where R > 2a. Equation (8) is interpreted as follows: the first term is just the rest mass of 
the P- or /-soliton and the remaining terms may be attributed to the interaction potential 
energy U{tq) of the impurity-soliton system, as a function of their separation ro- This attrac- 
tive interaction at short-range suggests that, in general, soliton solutions near nonmagnetic 
impurities may not be static. Spinless atoms may exert some force on the magnetic solitons, 
in addition to distorting them, and cause them to accelerate. If the force obtained from 
a potential agrees with the mass times acceleration of the soliton, then, due to the form 
of the effective potential in Eq.(8), which is an attractive well, we expect that the soliton 
center may have an oscillatory motion on the impurity. This Newtonian behavior may really 
happen to magnetic nonlinear excitations. For example, in 2D easy-plane magnetic systems, 
Wysin[^], [10| has shown that vortex excitations have motions that are the same as those 
predicted from a simple Newtonian dynamical equation of motion for the vortex center, 
involving a force F and effective mass M, i.e., 

? = M%, (10) 

where v is the velocity of the vortex center. In this case, the mass depends on the system 
size since a vortex is not a localized object. But, if the anisotropy is weak, the mass of the 



out-of-plane AFM vortices may actually reach a finite limit for large system size. Besides, 
it was shown that for a simple harmonic motion, the mass was found as 

M=~ 2 (11) 
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where uq is the translation mode frequency and K is the spring constant || 10| . Here, 



we assume that the results of Ref.@ extrapolate, with appropriate modifications, to the 
isotropic case, i.e., to the Lorentz- invariant 0(3) nonlinear a model, which represents the 
isotropic AFM. Thus, since the soliton is a localized structure, the rest mass would be 
expressed as M a = E a (R) / c 2 [[ll]], where c = 2JSa/h is the velocity of the long- wavelength 
spin waves (a =P or I). 

As it can be seen from Eq. (8), a small displacement tq of the P- or /-soliton center away 
from the exact impurity center (which is, in our model, the center of the circle of radius 
a) involves an energy increase. It is clear that these off-centered solutions are true bound 
states, because their energy is smaller than the free soliton mass. Any perturbation of the 
soliton will make it drift to the center of the well. The off-centered solutions are unstable. 
In fact, for displacements much less than the lattice spacing, the short-range potential well 
is close to harmonic, and can be described by an effective linear force acting on the soliton. 
The effective force constant K = (d 2 E /drQ) ro= o can be approximated by Eq. (9). Then, the 
frequencies of the oscillations of the solitons in a trapped state may be estimated analytically. 
Since this elastic constant is a function of E a (R), then, to determine the frequencies of the 
solitons oscillation modes, an analytical expression for the rest energy of the pinned-solitons 
would be useful. In Ref . || , it was shown that the energies of the two solitons can be written 
through the same expression E a (R) = E S R 2 /(R 2 + a 2 ), but with different a a 's (a p = 0.23a 
for the P and aj = 1.01a for the / soliton). Using these results, the solitons oscillation mode 
frequencies on a nonmagnetic impurity can be estimated as 



a 



2 / 2 2 \ 1 V 2 



^ 2 +(R 2 + a 



(12) 



_R 2 v a> Vv / 5[/? 4 -(2a) 4 ] 4(/? 2 

Figure 1 shows the frequencies for the two solitons as a function of their sizes. Of course, 
large solitons oscillate slower than small solitons. As a consequence, in general, P-solitons 
oscillate slower than /-solitons. 

The analysis of the spectrum of small-disturbance (the mesons of the theory or magnons) 
on a non- linear excitation always exhibits zero-frequency modes uj = 0, whenever the theory 



possesses a continuum symmetry. Although the soliton's energies (6) (or (8)) are invariant 
under spatial rotations around the z-axis (which for the hedgehog is equivalent to an iso- 
rotation around the internal 3-axis), they are not invariant under a translation in the x-y 
plane. In the presence of a nonmagnetic impurity, there is no possibility of a soliton to 
move away from a position by translations and keeping E unchanged. The impurity breaks 
the translational symmetry of the system, and thus, there are not zero-frequency modes 
associated to translations on the xy-plane (translation zero-modes). The oscillatory modes 
given by Eq.(12) can be considered to be driven by the force caused by the lattice defect (the 
spin vacancy), which can be evaluated independently of the spin- wave spectrum. To see this, 
we consider the change of the out-of-plane component of the Neel vector n% as the soliton 
center moves for a small displacement r << a from the impurity center. Since the soliton 
oscillates on the nonmagnetic impurity, we write r = Asin(u a t), where A is the oscillation 
amplitude assumed to be small. Then, the time dependent difference e a (r, t) = n 3a (r ) — 
^3«(0) gives the change of the out-of-plane component of the Neel vector as the soliton 
(P or i) oscillates. This difference can be easily evaluated. After linearizing in the small 
quantity, we obtain e a (r, t) ImlA^d^nsa exp(— iu a t)\. This expression shows that, if u a 
were zero, the function n 3a (A tl ) = n 3a (0) + A^d^n 3oi would describe a soliton which had been 
uniformly translated on the plane by a small amount A M = (A 1 , A 2 ) (translational mode). 
However, our solution has a nonzero frequency cu a , which implies in n 3a (A^) = n 3o (0) + 
A^d^n^a exp(— iuj a t) leading to an oscillatory motion for the soliton. These results can have 
a different interpretation. Here, the translational modes can be identified by viewing the 
motions of the spins that result when a given spin- wave eigenfunction is added to the original 
soliton structure. Then, the small-disturbance waveform e a (r, t) = A^d^n 3a exp(— icu a t) may 
also be viewed as out-of-plane spin waves with well defined frequency and angular momentum 
about the static soliton structure n 3a , which the center is localized exactly at the impurity 
center. It means that, instead of thinking about an oscillating soliton, we can imagine a 
stopped soliton surrounded by these magnons. In fact, using polar coordinates (r, 0) and 
supposing that A^ = (0, A), we have 



where m = ±1 is the angular momentum channel. In this interpretation, the above expres- 
sion gives the time dependent deformation of the static soliton. In particular, if one could 




(13) 
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observe the spin motions around the impurity (or in the central core of the soliton), their 
instantaneous directions could be used to estimate the position of the soliton center. Equa- 
tion (13) shows that the small-disturbance amplitude goes to zero as r — > oo and hence, it 
describes magnon modes, which are localized at a soliton. The soliton-impurity interaction 
may induce magnon local modes around the spin vacancy. These magnon local modes on 
the nonmagnetic impurity have frequencies that depend on the soliton type (P or / ) and 
size. In a discrete lattice, the possible sizes R are given by Ri = la (I is a positive integer) 
and then, the possible modes have also a discrete spectrum, whose frequencies are 



u a ,i = y/2- 
a 



S ~ /,2 ^ 1 1 

+ (I 2 + % 



1/2 



(14) 



I 2 " aJ VV5(Z 4 -16) 4(/ 2 + l) 2 , 

where 5 a is a number (5 P = 0.23 for the P-soliton and Si = 1.01 for the /-soliton). Note 
that I > 3, since the vacancy at the soliton center forbids the formation of solitons with size 
smaller than 3a in a discrete lattice. For I = 3, Up^ = 0.302c/a and ouj^ = 0.561c/a, while 
for I = 4, co>p )4 = 0.197c/a and uj^ = 0.402c/a, which typically have values (10 11 — 10 13 )s _1 . 
It is expected that local modes with well defined frequencies are induced around nonmag- 
netic impurities via soliton. Resonances at the characteristic frequencies can be in principle 
observed in electron spin resonance or inelastic neutron scattering experiments. Since these 
modes can occur only if P- or /- solitons are centered on the nonmagnetic impurity, reso- 
nance and neutron scattering measurements provide some methods to experimentally detect 
this effect and consequently, these two pinned solitons. 

Equation (12) keeps certain similarities with the results obtained by Wysin|9|, [H| for 



vortices in easy-plane magnets. The spin configuration of a vortex (or the magnetization 
inhomogeneity) extends to the limits of the system, causing its inertia against motion inside 
the system to depend on the system size. Hence, the effective vortex mass increases as 
the square of the system radius and the translational mode frequencies diminish as the 
reciprocal system size||, |10| . In our case, since the magnetization inhomogeneity depends 
on the soliton radius, the translational frequencies diminish as the reciprocal soliton size. 

The above results are valid only for small displacements satisfying ro << a. For large 
displacements, in which the soliton has to move for more than one lattice spacing, ro > a, 
nonlinear terms are present in the interaction force and besides, in real magnetic materials, 
discreteness effects must be considered. These effects can be described through a periodic 
pinning potential. These calculations are out of the scope of this paper. 
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Layered magnetic compounds of the form (C n i?2n+ 1-^-^3)2 MX4, where M is a transition 
metal ion and X is CI or Br, are nearly-classical spin systems for M = Mn 2+ and X — CI 
(the Mn ion has a spin 5/2). Hence, under nonmagnectic doping ions (such as Mg or Cd), 
these materials are good candidates to experimentally detect these modes and testing the 
theory. Finally, we would like to point out that our theory may also have some relevance 
to quantum spin systems [T2"|, |T3| , H |. In fact, the well known technique of spin coherent 



states effectively replaces spin operators by classical vectors S = (sin 6 cos 4>, sin 6 sin 0, cos 6) 
and incorporates the quantum features by means of the path integral over all space-time 
configurations of (6, <fi). Important 2D systems exhibiting quantum magnetism are the high- 
T c cuprate superconductors. In particular, concerning the problem of impurities, recent 
nuclear magnetic resonance measurements involving the high-T c cuprate superconductors, 
have shown that a substitution of the spin-1/2 Cu 2+ in the Cu — O plane by a strong 
nonmagnetic impurity, such as spin-0 Zn 2+ , enhances the antiferromagnetic correlations 



around the impurity [13]. We also suggest that these results have some relevance in the 



quantum Hall effect, where spin textures called skyrmions play important roles in the two- 



dimensional electron gas [14, 15 
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FIG. 1: Oscillation frequencies as a function of the soliton size. Solid curve is for P-solitons and 
dashed curve is for I-solitons. 



11 



